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ABSTRACT 



Minimal supersymmetric SU(5) GUTs are being squeezed by the recent values of a s , sm 2 9w, 
the lower limit on the lifetime for p —>■ vK decay, and other experimental data. We show 
how the minimal flipped SU(5) GUT survives these perils, accommodating the experimental 
values of ot s and sm 2 0w and other constraints, while yielding ap-* e/fi + n° lifetime beyond 
the present experimental limit but potentially accessible to a further round of experiments. We 
exemplify our analysis using a set of benchmark supersymmetric scenarios proposed recently 
in a constrained MSSM framework. 
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One of the key pieces of circumstantial evidence in favour of grand unification has long 
been the consistency of the gauge couplings measured at low energies with a common value 
at some very high energy scale, once renormalization effects are taken into account. This 
consistency is significantly improved when light supersymmetric particles are included in the 
renormalization-group running, in which case the agreement improves to the per-mille level IQ . 

However, this circumstantial evidence is not universally accepted as convincing. For exam- 
ple, it has recently been suggested that the logarithmic unification of the gauge couplings is as 
fortuitous as the apparent similarity in the sizes of the sun and moon 0. Alternatively, it has 
been argued that the unification scale could be as low as 1 TeV, either as a result of power-law 
running of the effective gauge couplings in theories with more than four dimensions ||, or in 
theories with many copies of the SU(3) x SU(2) x U{1) gauge group in four dimensions [§]. 

For some time now, detailed calculations have served to emphasize || how much fine tuning 
is needed in models with power-law running to reproduce the effortless success of supersym- 
metric grand unification with logarithmic running of the gauge couplings. Moreover, data from 
particle physics and cosmology provide independent hints for low-energy supersymmetry. Pre- 
cision electroweak data favour quite strongly a low-mass Higgs boson ||, as required in the 
minimal supersymmetric extension of the Standard Model (MSSM) f7|, and the lightest super- 
symmetric particle is a perfect candidate || for the cold dark matter thought by astrophysicists 
to infest the Universe. Many studies have shown that these and other low-energy data - such 
those on b — > 57 decay || and — 2 ||10| - are completely consistent with low-energy super- 



symmetry, and a number of benchmark supersymmetric scenarios have been proposed [i I 



Issues arise, however, when one considers specific supersymmetric grand unified theories. 
One is the exact value of sin 2 $w, which acquires important corrections from threshold effects 



at the electroweak scale, associated with the spectrum of MSSM particles [12], [13], and at the 



grand unification scale, associated with the spectrum of GUT supermultiplets |12|, [Mj. Precision 
measurements indicate a small deviation of sin 2 6y/ even from the value predicted in a minimal 
supersymmetric SU(5) GUT, assuming the range of a s [Mz) now indicated by experiment fl5 |. 



The second issue is the lifetime of the proton. Minimal supersymmetric SU(5) avoids the 
catastrophically rapid p — > e + 7r° decay that scuppered non-supersymmetric SU(5). However, 
supersymmetric SU(5) predicts p — > vK + decay through d = 5 operators at a rate that may be 
too fast [IE] to satisfy the presently available lower limit on the lifetime for this decay |T7], 



The latter requires the SU(5) colour-triplet Higgs particles to weigh > 7.6 x 10 16 GeV, whereas 
conventional SU(5) unification for a s (M z ) = 0.1185 ± 0.002, sin 2 6 w = 0.23117 ± 0.00016 and 
OLem{M z ) = 1/(127.943 ± 0.027) [0] would impose the upper limit of 3.6 x 10 15 GeV at the 
90% confidence level ||16|1 . This problem becomes particularly acute if the sparticle spectrum is 



relatively light, as would be indicated if the present experimental and theoretical central values 
of — 2 [[K|] remain unchanged as the errors are reduced. 



The simplest way to avoid these potential pitfalls is to flip SU(5) |L9|, P0|| . As is well 
known, flipped SU(5) offers the possibility of decoupling somewhat the scales at which the 
Standard Model SU(3), SU{2) and U(l) factors are unified. This would allow the strength 
of the U(l) gauge to become smaller than in minimal supersymmetric SU(5), for the same 



1 



value of a s (Mz) |13|]. Moreover, in addition to having a longer p — > e//i + 7r° lifetime than non- 
supersymmetric SU(5), flipped 577(5) also suppresses the d — h operators that are dangerous 
in minimal supersymmetric SU(5), by virtue of its economical missing-partner mechanism fT9| . 



In this paper, we re-analyze the issues of sin 2 6w and proton decay in flipped 577 (5) |13j , 
in view of the most recent precise measurements of a s (M z ) and sm 2 9w, and the latest lim- 
its on supersymmetric particles. We study these issues in the MSSM, constraining the soft 
supersymmetry-breaking gaugino masses m\/2 and scalar masses m to be universal at the 
GUT scale (CMSSM), making both a general analysis in the (mi/2,^0) plane and also more 
detailed specific analyses of benchmark CMSSM parameter choices that respect all the available 
experimental constraints [|TTJ. We find that the p — > e//i + 7r° decay lifetime exceeds the present 



experimental lower limit |L7j , with a significant likelihood that it may be accessible to the next 
round of experiments [^T|. We recall the ambiguities and characteristic ratios of proton decay 
modes in flipped SU(5). 

We first recall the lowest-order expression for a s (M z ) in conventional SU(5) GUTs, namely 

7 -a 

a s (M z ) = —h 7 • (1) 

5 sm u\y — 1 

The present central experimental value of a s {Mz) = 0.118 is obtained if one takes sin 2 0\y — 
0.231 and a -1 = 128, indicating the supersymmetric grand unification is in the right ball-park. 
However, at the next order, one should include two-loop corrections 52ioop as well as electroweak 
and GUT threshold corrections, that we denote by S\ ig ht and <5heavy Their effects can be included 
by making the following substitution in ([]]) [ |i~2]| : 

Sin 2 W — > Sin 2 9\y — f^loop — ^Hght — ^heavy , (2) 

where 5 2 ioo P ~ 0.0030, whereas <5i ig ht and 5heav y can have either sign. If one neglects 5n g ht 
and 5heavy, the conventional SU(5) prediction increases to a s (M z ) ~ 0.130 [0. A value of 



a s (M z ) within one standard deviation of the present central experimental value requires <5i ig ht 
and/or <5h ca vy to be non- negligible, so that the combination (52i oop + flight + ^heavy) is suppressed. 
However, in large regions of parameter space <5ii g ht > 0, which does not help. Moreover, in 
conventional SU(5), as was pointed out in |T2], 0, a compensatory value of 5h ea vy is difficult to 
reconcile with proton decay constraints. This problem is exacerbated by the most recent lower 
limit on r(p -> i)K + ) (T7| Q. 



As has been advertized previously fI5|| , an alternative way to lower a s (Mz) is to flip SU(5). 
In a flipped SU(5) model, there is a first unification scale M 32 at which the SU(3) and SU{2) 
gauge couplings become equal, which is given to lowest order by 

= *ib,*i, (3) 
2vr M z ' y J 
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. ^1„^, (4) 
2tt M z ' V ; 

1 It is true, as pointed out recently [Q, that this is not a problem if one allows arbitrary squark mixing 
patterns. However, such options must respect low-energy flavour-changing neutral-interaction limits |p3| , and 
are not possible in the CMSSM. 
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where a 2 = a/ sin 9wi a 3 — o s (M z ), and the one-loop beta function coefficients are b 2 = +1, 
6 3 = —3. The hypercharge gauge coupling ay = cos 2 9w) has, in general, a lower value a[ 
at the scale M 32 : 

a Y a[ 2tt M z ' w 

where by = 33/5. Above the scale M 32 , the gauge group is the full SU(5) x U(l), with the 
£7(1) gauge coupling a± related to a[ and the SU(5) gauge coupling a$ as follows: 

^ = l + ?i. (6) 

a x as ai 

The SU(5) and U{1) gauge couplings then become equal at some higher scale M51. The 
maximum possible value of M 32 , namely M.^ ax , is obtained by substituting a[ = «5(M 32 ) into 
(§D, and coincides with the unification scale in conventional SU(5): M^ ax = x exp((3 — 
8sin 2 ^vi/)vr/14a em (M z )), where M z = 91.1882 ± 0.0022 GeV |EJ. In general, one has 

as{Mz) = 5sin 2 ^-l + |aln(M 3 T x /M3 2 )' (?) 

and the flipped SU(5) prediction for a s (M z ) is in general smaller than in minimal SU(5), for 
the same value of sin 2 9w- The next-to-leading order corrections to (|7|) are also obtained by the 
substitution in @. Numerically, an increase of ~ 10% in the denominator in ([]]), which would 
compensate for the decrease due to <5 2 i 00 p, could be achieved simply by setting M 32 ~ |M.g ax 
in®. 

In order to understand the implications for r{p — > e//i + 7r°) decay, we first calculate M 32 , 
using (0) with sin 2 9w replaced by sin 2 6^ — #2ioop, leaving for later discussions of the possible 
effects of 5ii g ht,heavy- Fig- ID exhibits the correlation between M 32 and a s {M z ) in flipped SU(5). 
The solid lines indicate the range of values of M 32 allowed for a given value of a s (M z ) (as given 
in the MS prescription), assuming the experimentally-allowed range sin 2 9^ s = 0.23117 ± 
0.00016 ||18|| , and making no allowance for either light or heavy thresholds. For the central 



experimental value a s (M z ) = 0.1185, we see immediately that M 32 is significantly lower than 
its maximum value, which is M^ ax = 20.3 x 10 15 GeV for our central values of a s (M z ) and 
sin 2 9 W . 



We now explore the possible consequences of <5i ig ht for M 32 , following |T2], |13[. We approxi- 
mate the 5ii g ht correction by 



a 



20tt 



28 _ , 32 



flight = 7^ 3L(m t ) + —L(rrig) - — L(m^) - L(m h ) - AL(m H ) 



3 v y/ 3 



5 35 19 

-L(m 9 ~) - 3L(m iL ) + 2L(m ia ) - ^ L ( m t 2 ) - ^ 



where L(x) = ln(x/M z ). As already mentioned, we assume that the soft supersymmetry- 
breaking scalar masses m , gaugino masses my<z and trilinear coefficients A are universal at 
the GUT scale (CMSSM). We used ISASUGRA |25j to calculate the sparticle spectra in terms of 



3 




Figure 1: The solid lines show the correlation between M 32 in flipped SU(5) and a s {Mz) in the 
MS prescription, assuming sin 2 9yf s = 0.23117 ± 0.00016, including 52ioop but neglecting <5i ig ht 
and 5heavy The points indicate the changes in r(p — > e/ fi + n°) found for a s {Mz) = 0.1185 and 
the central value of sin 2 9w when including also the values of 5u g ht calculated for the CMSSM 
benchmark points. 

these quantities, tan/3 and the sign of fi, assuming m t = 175 GeVQ In evaluating (g), (mn) 
{rriq) {mi) were interpreted as the geometric means of the chargino and neutralino (H, A, H^) 
(m, d, s, c) (e, jl) masses, respectively, and the mixings of f, b and t were all taken separately 
into account. 

The unknown parameters in (|8]) were constrained by requiring that electroweak symmetry 
breaking be triggered by radiative corrections, so that the correct overall electroweak scale 
and the ratio tan/3 of Higgs v.e.v.'s fix and in terms of mi/2 and rriQ. Before making 
a more general survey, we recall that a number of benchmark CMSSM scenarios have been 
proposed |TTJ, which include these constraints and are consistent with all the experimental 



limits on sparticle masses, the LEP lower limit on m/j, the world-average value of b — > S'j 
decay, the preferred range 0.1 < Q x h 2 < 0.3 of the supersymmetric relic density, and g M — 2 
within 2 o~ of the present experimental value. These points all have A Q = 0, but otherwise 
span the possible ranges of m^, m , tan/3 and feature both signs for /i. Fig. [I] also shows the 
change in M32 induced by the values of 5n g ht in these benchmark models, assuming a fixed 
value a s (Mz) = 0.1185. In general, these benchmark models increase M32 for any fixed value 
of a s (Mz) and sin 2 #w As a s (Mz) varies, the predicted value of M 32 in each model varies in 
the same way as indicated by the sloping lines. We recall that the estimated error in a s [Mz) 
is about 0.002, corresponding to an uncertainty in M 32 of the order of 20%, and hence a 



2 Heavy singlet neutrinos were not used in the renormalization-group equations. 
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corresponding uncertainty in the proton lifetime of a factor of about two. The error associated 
with the uncertainty in sin 2 Ow is somewhat smaller 

We now turn to the calculation of r(p — ► e/[i + n°). We recall first that the form of the 
effective dimension-6 operator in flipped SU(5) is different |24 , |26| from that in conventional 
^(5)0,0: 

9l 



e^dle^^PLd^ua^L) + h.c. 

32 

ijkfjc 2irju n , -c 2ir/ 2 i 



+ (e» k (d c k e 2 ^ cos# c + s c k e^ sin e e )-fP LUj )(uavP L £ L ) + h - c ] ( 9 ) 

where 9 C is the Cabibbo angle Q. Also appearing in @ are two unknown but irrelevant CP- 
violating phases 7711,21 and lepton flavour eigenstates and II that are related to mass eigen- 
states by unknown but relevant mixing matrices: 

v L = v F U v , £ L = £ F U t . (10) 
Despite our ignorance of the mixing matrices (|T0|), some characteristic flipped SU (5) predictions 



can be made 24 



r(p - e + TT ) = _^|C/ €ii |2 r(p ^ p7r+) = cos 2^|^ n |2 r(n ^ ^ 

r(n e+7r~) = 2T(p -> e + 7r°) , r(n -> // + 7r") = 2r(p -> /i+7r°) 
r(p - = ^-^\U il2 \ 2 T(p - Z/7T+) = cos 2 £ c |f/, 12 | 2 i> - z7vr°) (11) 

In the light of recent experimental evidence for near-maximal neutrino mixing, it is reasonable to 
think that (at least some of) the e/fi entries in Ui are 0(1). In what follows, we assume that the 
lepton mixing factors |t^ 1112 | 2 are indeed 0(1), and do not lead to large numerical suppressions 
of both the p — > e//i + 7r° decay rates. Note that there is no corresponding suppression of the 
p — ► utt + and n — > z/7r° decay rates, since all the neutrino flavours are summed over. However, 
without further information, we are unable to predict the ratio of p — > e + X and p — > fi + X 
decay rates. Hereafter, wherever we refer to p — > e + 7r° decay, this mixing-angle ambiguity 
should be understood. 

The p — > e + 7r° decay amplitude is proportional to the overall normalization of the proton 
wave function at the origin. The relevant matrix elements are a, (5, defined by 

(0|e iifc (^') fl «||p(k)) = aur(k), (12) 
(0|e^(^) L <|p(k)) ee /3 Ui (k). (13) 



3 We note from Fig. [l] that there is no benchmark model for which conventional SU (5) grand unification is 
possible, with the measured values of a s (Mz) and sin 2 ^^; unless one invokes GUT threshold effects. 

4 Note the absence [^J], ^(| in the corresponding decay rate of the factor (1 + (1 + V^l 2 ) 2 ) found ]27], |2^| in 
conventional SU(5), as recently re-emphasized in [[t6|. This lengthens t p by ss 5 in flipped SU(5), an effect that 
is typically more than offset by the reduction in M32. 
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The reduced matrix elements a, {3 have recently been re-evaluated in a lattice approach p9 
yielding values that are very similar and somewhat larger than had often been assumed previ- 
ously, and therefore exacerbating the proton-stability problem for conventional supersymmetric 
SU(5). Here, we use here the new central value a = (3 = 0.015 GeV 3 for reference. The error 
quoted on this determination is below 10%, corresponding to an uncertainty of less than 20% in 
r(p — > e + 7r°), which would be negligible compared with other uncertainties in our calculation. 
Thus, we have the following estimate, based on |26|, [16| and references therein: 




r(p- e V) = 3.8 x 10* — R^f^ 7 ( 14 ) 



for use in the subsequent analysis, where we have absorbed reference values for M 32 and a^(M 32 ) 
as well as a and p, and a 5 (M™ x ) / a 5 (M 32 ) = 1 - (33/28) (a 5 (M™ x )/2ir)\n(M 32 /M. 



max") 
32 ) 



We present a general view of flipped SU(5) proton decay in the CMSSM in Fig. |2|. The 
thick solid (blue) lines are contours of r(p — > e + 7r°) for the indicated choices of tan/5 and 
the sign of \x Q, which span (most of) the range of possibilities. Where applicable, we have 
indicated by (blue) crosses and labels the CMSSM benchmark points with the corresponding 
value of tan/3 and sign of ji [[11]]. Following |]30| , the dark (red) shaded regions in the bottom 



right-hand parts of each panel are excluded because the LSP is the lighter f: astrophysics 
excludes a charged LSP. The light (turquoise) shaded regions have LSP relic densities in the 
preferred range 0.1 < Q x h 2 < 0.3 for cold dark matter. The intermediate (green) shaded 
regions at lower mm are excluded by b — > 57, which is a more important constraint for fi < 0. 
The other shaded (pink) regions at large (m^, mo) are consistent with — 2 at the 2-a 
level. In panels (c) and (d), the hatched regions at low my 2 and large m are those where 
electroweak symmetry breaking is no longer possible, and the horizontally-striped regions at 
low mo have tachyons. The dash-dotted (blue) line at small (mi/2,mo) i n panel (a) corresponds 
to rrie = 100 GeV. The near-vertical dashed (black) lines at small m.1/2 correspond to the LEP 
lower limit m^ = 103.5 GeV, and the dot-dashed (red) lines to LEP lower limit m^ = 114 GeV 



as calculated using the FeynHiggs code pifl . In each case, only larger values of m 1 / 2 are allowed, 
although there is uncertainty in the location of the m^ line []. 

We see in Fig. I that the 'bulk' regions of the parameter space preferred by astrophysics 
and cosmology, which occur at relatively small values of (m.1/2, mo), generally correspond to 
r(p — > e + 7r°) ~ (1 — 2) x 10 35 y. However, these 'bulk' regions are generally disfavoured by 
the experimental lower limit on rrih and/or by b — > sj decay. Larger values of r(p — > e + 7r°) are 
found in the 'tail' regions of the cosmo logical parameter space, which occur at large mi/2 where 
X~ £ coannihilation may be important, and at larger m.1/2 and mo where resonant direct-channel 
annihilation via the heavier Higgs bosons A, H may be important. 

We turn finally to the possible implications of the GUT threshold effect <5heav y [P| Ell - ^ 



5 The horizontal spacing between points sampled was comparable to the thickness of these lines. 

6 For fuller discussions of the implementations of these constraints with and without ISASUGRA, see jllL 
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Figure 2: TTie solid (blue) lines are contours of r(p — > e/// + 7r°) in the (mi/2,mo) plane for 
the CMSSM with (a) tan/? = 10,// > ; (6j tan/? = 10,// < 0, fcj tan/? = 35,// < 
and C<ij tan/3 = 50,// > 0. TTie (blue) crosses indicate the CMSSM benchmark points with the 
corresponding value of tan. /3 and sign of fx Following j^Jffl, the dark (red) shaded regions are 
excluded because the LSP is charged, the light (turquoise) shaded regions have 0.1 < Q x h 2 < 0.3, 
intermediate (green) shaded regions at low mu2 are excluded by b — > s'-f, shaded (pink) regions 
at large (my2, m o) ore consistent with g^ — 2 at the 2-a level, and electroweak symmetry breaking 
is not possible in the hatched regions. The near-vertical dashed (black) lines correspond to the 
LEP lower limit m£ = 103.5 GeV, the dot-dashed (red) lines to = 114 GeV as calculated 
using the FeynHiggs code [31] , and the dotted (blue) lines at small (jtoi/2, tuq) to = 100 GeV. 
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general expression for this in flipped SU(5) is given in 



12 



a 

' heavy 



20tt 



, M 32 , M 32 , M 32 
-6 In — — - 6 In — — + 4 In — p- 

M Hs M Bs My 



a 

20^ 



A4A5 



(15) 



where Mh 3 = X^\V\ and Mg 3 = As| V| are the masses of the heavy triplet Higgs supermul- 
tiplets, the X,Y gauge bosons and gauginos have common masses My = gs\V\ where V is 
the common v.e.v. of the 10 and 10 Higgs supermultiplets, A^s are (largely unconstrained) 
Yukawa couplings, g 5 is the 577(5) gauge coupling, and r = max{g 5 , A 4 , A 5 }. Thanks to the 
economical missing-partner mechanism of flipped 577(5), the H 3 and H 3 do not mix, and hence 
do not contribute significantly to proton decay. Thus there is no strong constraint on M Ha g 3 
from proton decay in flipped 577(5), and it is possible that M Ha g 3 < My (i.e., r = g 5 ). In 
this case, we can see from flTBD that ^heavy < naturally For instance, as pointed out in |13 



if A4, A5 ~ hg 5 , then <5heav y ~ —0.0030, which completely compensates the <5 2 i 00 p contribution. 



We also recall that, in general, including £heavy leads to a re-scaling of the M 32 /M: 



max. 
32 • 



Mi ax M$ ax 

We display in Fig. |3] the possible numerical effects of <5heavy on r{p — > e/yU + 7r°) in the various 
benchmark scenarios, assuming the plausible ranges —0.0016 < <5h ca vy < 0.0005 [[y|. The 
boundary between the different shadings for each strip corresponds to the case where 5h ea v y = 0. 
The left (red) parts of the strips show how much r(p — > e + n°) could be reduced by a judicious 
choice of <5heavy, and the right (blue) parts of the strips show how much r(p — > e + n°) could 
be increased. The inner bars correspond to the uncertainty in sin 2 0yy. On the optimistic 
side, we see that some models could yield r(p — > e + 7r°) < 10 35 y, and all models might have 
r(p — > e + 7T°) < 5 x 10 35 y. However, on the pessimistic side, in no model can we exclude the 
possibility that r(p — ► e + 7r°) > 10 36 y. 

We recall that a new generation of massive water-Cerenkov detectors weighing up to 10 6 tonnes 
is being proposed pi| , that may be sensitive to r(p — > e + 7r°) < 10 35 y. According to our calcu- 



lations, such an experiment has a chance of detecting proton decay in flipped 577(5), though 
nothing can of course be guaranteed. We recall that there is a mixing-angle ambiguity ( |TTD 
in the final-state charged lepton, so any such next-generation detector should be equipped to 
detect e + and/or ji + equally well. We also recall |24], |26| that flipped 577(5) makes predictions 
([TTD for ratios of decay rates involving strange particles, neutrinos and charged leptons that 
differ characteristically from those of conventional SU(5i). Comparing the rates for e + , /i + and 
neutrino modes would give novel insights into GUTs as well as mixing patterns. 

We conclude that flipped SU(5) evades two of the pitfalls of conventional supersymmetric 
577(5). As we have shown in this paper, it offers the possibility of lowering the prediction for 
a s [Mz) for any given value of sin 2 6w and choice of sparticle spectrum. As for proton decay, we 
first recall that flipped 577(5) suppresses p — > uK + decay naturally via its economical missing- 
partner mechanism. As in conventional supersymmetric SU{5), the lifetime for p — > e/fi + 7i° 
decay generally exceeds the present experimental lower limit. However, as we have shown in 
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1 10 35 2 10 35 5 10 35 110 36 2 10 36 5 10 36 
X P [Y] 

Figure 3: For each of the CMSSM benchmark points, this plot shows, by the lighter outer bars, 
the range ofr(p — > e//i + 7r°) attained by varying 5h cscvy over the range -0.0016 to + 0.0005 ]I3|/ . 
The central boundary of the narrow inner bars (red, blue) corresponds to the effect ofSng^ alone, 
with 5heavy = 0, while the narrow bars themselves represent uncertainty in sin 2 6 W . We see that 
heavy threshold effects could make r{jp — > e/p + 7t°) slightly shorter or considerably longer. 

this paper, the flipped 577(5) mechanism for reducing a s (Mz) reduces the scale M32 at which 
colour SU(3) and electroweak ST/ (2) are unified, bringing r(p — > e/fi + ir ) tantalizingly close 
to the prospective sensitivity of the next round of experiments. Proton decay has historically 
been an embarrassment for minimal SU(5) GUTs, first in their non-supersymmetric guise and 
more recently in their minimal supersymmetric version. The answer may be to flip SU(5) out 
of trouble. 
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